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PREFACE 

Authority  for  the  investigation  reported  herein  is  contained  in  FY  1961 
Instructions  and  Outline,  Military  Construction  Investigations,  Engineering 
Criteria  and  Investigations  and  Studies,  Investigation  of  Arctic  Construction, 
Creep  of  Frozen  Soils. 

This  study  was  conducted  for  the  Engineering  Division,  Directorate  of 
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ministered  by  the  Civil  Engineering  Branch,  Mr.  T.B.  Pringle,  Chief. 

Mr.  Francis  H.  Sayles,  Research  Civil  Engineer,  Applied  Research 
Branch  carried  out  the  study  and  prepared  this  report.  The  investigation 
was  under  the  general  direction  of  Mr.  K.A.  Linell,  Chief,  Experimental 
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Laboratory  (CRREL),  U.S.  Army  Terrestrial  Sciences  Center  (USA  TSC). 
Personnel  assisting  in  the  investigation  were  SP  Richard  Putnam  and  SP 
Richard  O.  Lunde.  Mr,  Robert  Bonnett  assisted  in  the  testing. 

This  report  has  been  critically  reviewed  by  Professor  Clyde  E.  Kesler 
of  the  University  of  Illinois  and  by  Mr.  Frederick  J.  Sanger  of  USA  TSC. 

The  author  wishes  to  thank  the  reviewers  and  Dr.  A.  Assur,  Chief  Scientist, 
USA  TSC,  for  iheir  constructive  suggestions. 

Lieutenant  Colonel  John  E.  Wagner  was  Commanding  Officer /Director  of 
the  U.S.  Army  Terrestrial  Sciences  Center  during  the  publication  of  this  re¬ 
port  and  Mr.  W.  K.  Boyd  was  Chief  Engineer. 

USA  TSC  is  a  research  activity  of  the  Army  Materiel  Command. 
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SUMMARY 


Unconf  ned  compressive  cieep  strengths  ana  strains  were  measured 
for  frozen  saturated  Ottawa  sand  (20-30)  and  Manchester  fine  sand.  The 
creep  tests  were  conducted  at  approximate  stress  levels  of  60,  35,  20  and 
5%  of  the  conventional  unconfined  compressive  strength.  Testing  tempera¬ 
tures  were  15,  25,  29  and  31F.  It  was  found  that  the  unconfined  compres¬ 
sive  creep  strength  of  the  frozen  sand  can  be  predicted  using  Vialov's 
strength  formula;  that  creep  strain  can  be  predicted  using  two  short-term, 

t4> 

high-stress-level  creep  tests  using  «  =  "ti  •  “  1  that  tota  strain  can  be 

1  1/K  t4' 

predicted  using  «  =  -  — —  +  t_  ;  and  that  for  stresses 

«?[«/«.)*  i]»J  *  " 


below  the  long-term  strength,  the  strain  rate  is  directly  proportional  to 
the  reciprocal  of  time  during  stress  action  until  complete  stabilization 
occurs,  (ij  =  strain  rate  1  nour  after  stress  is  applied;  t  =  time; 

4/  =  (M-l)/M,  where  M  =  o-1'w  and  w  is  a  constant  for  each  materia1; 


(r01  =  stress  at  9  =  0;  9  =  temperature  in  degrees  below  freezing  point  of 
water;  90  =  a  constant  reference  value  of  G;  a  and  K  are  constants; 

*0=  initial  instantaneous  strain. ) 
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INTRODUCTION 

The  design  of  stable  structures  on  permafrost  requires  a  knowledge  of 
the  strength  and  deformation  characteristics  of  frozen  soil.  Published  mat¬ 
erial  on  the  strength  and  deformation  properties  of  frozen  soil  prior  to  1952 
was  of  Russian  origin  and  was  generally  incomplete  as  to  description  of  soils 
and  testing  procedures.  In  1952  the  former  Arctic  Construction  and  Frost 
Effects  Laboratory  (ACFEL)*  of  the  U.  S.  Army  Engineer  Division,  New  Eng¬ 
land,  published  a  report  summarizing  experimental  data  obtained  up  to  that 
time,  including  the  results  of  their  investigations  (ACFEL,  1952).  Since  1952 
the  Russians,  notably  Tsytcvich  and  Vialov,  have  published  rather  complete 
experimental  data  on  the  strength  and  deformation  properties  of  some  naturally 
frozen  silts  and  clays  (Tsytovich,  1954,  1958;  Vialov,  1959;  Vialov  et  al. , 
1962;  Vialov  and  Tsytovich,  1955).  In  addition,  they  summarized  and  formu¬ 
lated  theories  and  empirical  equations  relating  strength  and  deformation  of 
frozen  soils  to  the  soil  temperature  and  duration  of  the  applied  load.  Sanger 
and  Kaplar  (1963)  published  deformation  data  and  empirical  equations  relat¬ 
ing  unconfined  compressive  deformation  and  rate  of  deformation  to  applied 
stress  and  temperature.  This  investigation  included  a  variety  of  Soils,  tested 
at  various  temperatures  from  about  18F  to  32F.  Each  creep  test  was  limited 
to  60  hours  duration. 

The  purpose,  of  this  investigation  is  to  evaluate  the  influence  of  tempera¬ 
ture  and  stress  on  creep  and  long-term  strength  of  saturated  frozen  sands, 
and  to  provide  data  for  design  in  frozen  soils. 

This  report  with  its  appendix  presents  the  completed  results  of  the  un¬ 
confined  compression  tests  performed  on  saturated  Ottawa  sand  (20-30)  and 
Manchester  fine  sand.  This  is  only  the  first  phase  of  the  current  investiga¬ 
tion  which  includes;  ( 1 )  unconfined  compression  creep  tests  on  Ottawa  sand, 
Manchester  fine  sand,  New  Hampshire  silt  and  a  clay  and  (2)  triaxial  creep 
testing  of  Ottawa  sand. 


DEFINITION  OF  TERMS 

Instantaneous  strength  is  the  maximum  stress  determined  by  loading  the 
test  specimen  at  a  constant  strain  rate  of  0.  033/min. 

Long-term  strength  is  the  maximum  stress  that  the  frozen  soil  can  with- 
stand  indefinitely  and  exhibit  either  a  zero  or  continuously  decreasing 
strain  rate  with  tirpe. 


*  ACFEL  was  merged  with  the  former  U.S.  Army  Snow,  Ice  and  Permafrost 
Research  Establishment  (SIPRE)  in  1961  to  form  the  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory. 
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Cun  vc. it  ion  ill  strain  is  the  axial  deformation  divided  by  the  original  length 
of  the  specimen  (ec  AL/L0  ). 

True  strain  is  the  axial  deformation  at  a  given  instant  of  time  divided  by  the 
actual  specimen  length  at  tha*  time  [in  terms  of  conventional  strain, 
e t  ln(  1  /  (1  • 

Incipient  failure  of  a  test  specimen  occurs  when  the  strain  rate  starts  to  in¬ 
crease  with  time  (start  of  tertiary  creep),  after  a  period  of  minimum 
strain  rate  (point  l.  Fig.  25). 

F ailure  of  a  compression  test  specimen  of  frozen  sand  means  a  continuous 

loss  of  resistance  to  loading,  after  reaching  a  peak.  It  occurs  by  either 
an  abrupt  brittle-type  fracture  or  a  plastic  flow  accompanied  by  fissur- 
ing  and  often  ending  in  rupture  of  the  specimen. 

Elastic  deformation  disappears  entirely  upon  release  of  Tie  stress  which 
caused  it. 


Delayed  elastic  deformation  is  elastic  deformation  that  requires  a  noticeable 
period  of  time  to  occur  or  recover.  Theoretically,  elastic  deformation 
velocities  slower  than  the  speed  of  spurd  are  delayed.  In  this  report 
delayed  elastic  recovery  is  slower  than  the  speed  of  sound  periods  by 
hours  and  days . 

Viscous  deformation  is  an  irreversible  deformation  in  which  the  rate  of  de¬ 
formation  depends  upon  the  applied  stress. 

Plastic  deformation  is  an  irreversible  deformation  which  is  independent  of 
time. 


Stress  ratio,  p 


Applied  constant  stress 
''instantaneous"  strength 


REVIEW  OF  THEORY 


Deformation 

Vialov  and  Tsytovich  (1955)  explain  the  physical  process  of  creep  in  frozen 
soil  by  considering  the  condition  of  applying  a  constant  lead  to  a  frozen  soil 
mass.  This  load  concentrates  the  stress  between  the  soil  particles  at  their 
points  of  contact  with  the  ice,  causing  pressure-melting  of  the  ice.  Differences 
in  water  surface  tensions  are  produced  and  the  unfrozen  water  moves  to  re¬ 
gions  of  lower  stress  where  it  refreezes.  The  process  of  ice  melting  and  water 
movement  is  accompanied  by  a  breakdown  of  the  ice  and  structural  bonds  of  the 
soil  grains,  the  plastic  deformation  of  the  pore  ice  and  a  readjustment  in  the 
particle  arrangement,  the  result  of  which  is  the  time-dependent  deformation 
phenomenon  of  creep.  This  structural  deformation  leads  to  a  denser  packing 
of  the  soil  particles,  which  in  turn  causes  a  strengthening  of  the  material  due 
to  the  increased  number  of  firm  contacts  between  soil  grains  and  hence  an  in¬ 
crease  ir  internal  friction  between  grains  (ACFEL,  1952).  During  this  process 
there  is  also  a  weakening  of  the  structural  cohesion  and  possibly  an  in¬ 
crease  in  the  amount  of  unfrozen  water  in  the  frozen  soil  (particularly  in  fine¬ 
grained  soils).  All  of  this  action  is  tirne-dependent.  If  the  applied  load 
does  not  exceed  the  long-tetm  strength  of  the  frozen  soil,  then  the  weakening 
process  is  compensated  by  the  strengthening;  the  deformation  is  damped,  i.e., 
the  rate  of  deformation  decreases  with  tim  -.  However,  if  the  applied  load  ex¬ 
ceeds  the  frozen  soil  long-term  strength,  the  breakdown  of  internal  bonds  is 
not  completely  compensated  by  the  strengthening  process  and  then  the  rate  of 
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deformation  increases  with  timei  resulting  in  undamped  deformation  which 
eventually  develops  into  plastic  flow  and  ends  in  a  breakdown  of  the  frozen 
soil  structure. 

Experiments  *  show  that  the  deformation  characteristics  of  frozen  soil 
are  similar  to  those  depicted  by  the  classical  creep  curve  for  metals  (Fig. 
1).  As  an  aid  in  studying  these  deformation  characteristics,  Vialov  and 
others  h.  ve  proposed  mechanical  rheological  models.  In  Figure  1,  Vialov's 
model  (1959)  with  the  various  components  labeled,  depicts  the  first  and 
second  stages  of  creep,  but  does  not  include  the  initial  plastic  deforma¬ 
tion  or  third  stage  of  creep,  i.  e.  ,  the  visco-plastic  flow  preceding  complete 
collapse  of  the  material  structure. 


RUPTURE 


o  limit 


MECHANICAL  MOOEL  A  MODIFIED  VIALOV  MECHANICAL  MODEL 

■kfttr  Vmlov 


Figure  1.  Mechanical  rheological  models. 


+  Vialov,  1962,  Fig.  27,  (Translation  page  lOi)  for  silt  and  clay  and  Figures 
14  and  15  of  this  report  for  saturated  Ottawa,  sand  and  Manchester  fine  sand. 
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In  Figure  1  Vialov's  model  is  modified  to  include  plastic  characteristics. 
The  various  elements  on  the  model  are  labeled  to  correspond  to  the  appro¬ 
priate  section  of  the  classical  creep  curve.  In  the  modified  model,  the  initial 
plastic  deformation  is  represented  by  a  frictional  slide.  An  air  gap  (this  gap 
may  have  finite  value  or  be  zero)  between  the  slide  and  the  Voigt  element  per¬ 
mits  displacement  to  take  place  before  viscoelastic  movement  begins.  If  the 
load  exceeds  the  yield  point  of  the  material,  the  slide  will  mo'-e  to  reduce  or 
close  the  gap,  thus  producing  instantaneous  plastic  deformation.  When  the 
load  is  removed,  this  part  of  the  deformation  is  not  recovered.  In  some  in¬ 
stances  elastic  deformation  may  occur  without  reaching  the  plastic  limit  upon 
loading,  thus  producing  purely  elastic  deformation.  The  viscoelastic  element 
is  the  Voigt  or  Kelvin  element  which  contributes  a  delayed  elastic  effect.  The 
dashpot  of  element  4  represents  the  viscous  portion  of  the  curve.  The  fric¬ 
tion  element  (blocking  device)  requires  the  activating  force  to  reach  a  certain 
magnitude  before  plastic  and  viscous  deformation  can  occur.  To  represent 
viscoplastic  flow  a  specially  shaped  dashpot  is  shown  to  permit  an  increasing 
rate  of  flow  as  the  material  approaches  failure.  The  modified  model  is  too 
complicated  for  a  practical  numerical  analysis  and  is  merely  shown  to  illus¬ 
trate  the  deformation  components. 


Strength  of  frozen  soils,  as  with  unfrozen  cohesive  soils,  depends  upon 
both  the  cohesion  and  the  internal  friction  of  the  component  materials.  In  fro¬ 
zen  soils  the  cohesion  component  according  to  Vialov  and  Tsytovich  (l°55) 
can  be  attributed  to:  (1)  the  molecular  forces  of  attraction  between  solid 
particles;  (2)  physical  or  chemical  cementing  of  particles  together;  ard 
(3)  cementing  the  soil  particle  by  ice  formation  in  the  soil  voids.  Cementing 
by  ice  is  the  result  of  the  bonds  between  the  ice  crystals  and  the  soil  particles 
even  though  the  soil  particles  aie  surrounded  by  a  film  of  unfrozen  'vater. 

This  unfrozen  water  is  under  the  influence  of  molecular  forces  of  the  soil 
particles  and  it  seems  possible  that  the  strongly  attached  water  molecules  are 
capable  of  transmitting  normal  and  shear  forces  between  solid  ice  and  solid 
grains.  The  ice  cohesion  depends  upon  the  amount  of  ice,  the  strength  of  the 
ice,  and  the  area  of  ice  in  contact  with  the  soil  particles,  each  of  which  de¬ 
pends  upon  the  soil  temperature.  The  internal  friction  depends  on  the  soil 
grain  arrangement,  sizes,  distribution,  shape,  and  on  the  number  of  grain-to- 
grain  contacts.  It  is  emphasized  that  ice  cohesion  is  the  dominant  strength 
factor  in  frozen  soil  even  though  internal  fricti  on  becomes  significant  in  dense 
sands. 


TESTING 


The  unconfined  compression  test  was  chosen  as  the  primary  test  for  this  inves¬ 
tigation  because  of  its  simplicity  and  its  suitability  for  adoption  as  a  field  laboratory 
test. 


In  addition  to  the  compression  test,  sonic  tests  and  ball  penetration  tests 
were  performed.  The  sonic  tests  were  exploratory  in  nature  with  the  primary 
purpose  of  testing  equipment  and  developing  techniques.  Both  tests  showed 
some  promise  but  the  allotted  time  did  not  permit  continuance  until  techniques 
could  be  perfected. 
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APPARATUS 


Freezing  facilities 

Soil  and  ice  specimens  were  frozen  in  freezing  cabinets  mounted  in  a 
USA  CRREL  walk-in  tvpe  cold  room  maintained  at  +  40F  +  1F.  The  freez 
ing  cabinets  used  in  this  project  (Fig.  2,  3)  were  equipped  with  hinged  covers 
on  top  and  a  thermal-pane  window  in  the  front.  Insulation  was  provided  in 
the  sides  and  cover.  The  bottom  of  the  cabinet  consisted  of  at.  expanded 
metal  grill  allowing  the  bottom  of  the  specimens  to  be  exposed  to  the  40F 
-oom  temperature  during  freezing  while  the  tops  of  the  specimens  were 
subjected  to  the  desired  freezing  temperature.  The  freezing  cabinet  was 
cooled  by  coils  mounted  in  its  sides  and  back  wall,  from  13  in.  above  the 
bottom  to  the  top  of  the  cabinets. 

The  cabinet  temperature  could  be  controlled  to  wituin  ±  0.  5F  by  means 
of  heating  coils  located  in  the  air  stream  of  the  air  circulation  fan.  The 
heat  supplied  was  regulated  by  a  Bayley  temperature  controller  Model  122. 
De-aired  water  was  supplied  to  the  bottom  of  the  specimens  during  freezing 
by  means  of  an  external  reservoir. 


Figure  2.  View  of  four  freezing  cabinets  inside 
40F  cold  room. 
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Figure  3.  Front  of  freezing  cabinet  showing  freezing  mold 
in  position  with  thermocouple  leads  to  the  center  specimen. 
The  bimetallic  thermal-regtflator  is  shown  to  right  of  view¬ 
ing  window. 


F reezing  mold 

The  freezing  mold  (Fig.  4,  5)  consisted  of  a  Plexiglas  block  17  1/4  in. 
square  and  7  in.  thick,  through  which  25  3-in.  diam  holes  were  bored.  Each 
hole  was  fitted  with  a  split  sleeve  having  a  wall  thickness  of  about  1/8  in. 
to  permit  specimen  ejection  from  the  mold  without  subjecting  the  specimen 
to  the  ejection  force.  The  top  and  bottom  of  the  mold  were  covered  by  1/2-in. 
thick  aluminum  plates  sealed  to  the  meld  block  with  1 /2-in.  soft  rubber  gas¬ 
kets.  When  the  freezing  mold  was  assembled,  an  expanded  metal  screen  with 
1/2x1  /4-in.  openings,  a  200-mesh  bronze  screen  and  a  muslin  mat  were 
placed  at  the  top  and  bottom  of  the  mold  block  to  retain  the  unfrozen  soil  speci¬ 
mens  in  the  cylinders  and  to  act  as  a  filter. 

Loading  equipment 

Three  types  of  loading  devices  were  used  in  the  testing  program  to  accom¬ 
modate  the  different  strength  and  deformation  characteristics  of  the  frozen 
sands. 

Tests  to  determine  the  instantaneous  compressive  strength  and  short-term 
creep  strength  of  frozen  soils  with  relatively  high  resistance  were  performed 
in  a  20,  000-lb  capacity  air-actua^d  hydraulic  press  (Fig.  6).  Loads  were 
applied  to  the  test  specimen  by  means  of  an  unconfined  test  chamber  placed  in 
the  press.  The  press  is  capable  of  head  movement  rates  up  to  18  in.  /min. 

For  creep  tests,  vibration-free  constant  loads  can  be  maintained  by  the  hy¬ 
draulic  press  for  extended  periods  ot  time  within  2%  of  the  applied  load  for 
loads  greater  than  3000  lb. 
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Figure  4a.  Empty  freezing  mold  showing  projecting  split  sleeves 


! 


Figure  4b.  Freezing  mold  filled  with 
sleeves.  Mold  is  supported  on 


dry  sand  showing  split 
ejection  frame. 


Figure  6.  Pneumatically  actuated  hydraulic  press. 


CREEP  OF  FROZEN  SANDS 


Figure  5.  Freezing  mold  charged  with  specimens  being  de- 
aired  and  saturated. 
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Figure  7.  Constant  stress  appa¬ 
ratus  showing  programming  cam 
and  load  measuring  system. 


Creep  tests  in  which  large  deforma¬ 
tion  occurs  were  performed  on  a  con¬ 
stant  stress  press  (capacity  4000  lb)(Fig. 

7) .  This  press  features  a  programming 
cam  that  maintains  constant  axial  stress 
to  within  1%  of  applied  stress  on  the  test 
specimen  during  deformation.  The  load¬ 
programming  is  based  on  the  assump¬ 
tions  that  the  cross-sectional  area  of  the 
specimen  remains  uniform  throughout  its 
length  during  deformation,  and  that  the 
volume  remains  constant  throughout  the 
test. 

Long-term  creep  tests  resulting  in 
small  deformation  were  performed  on  a 
lever-type  press  (capacity  2000  lb)  (Fig. 

8) .  As  the  sample  deforms,  the  height 
of  the  fulcrum  is  adjusted  to  maintain  the 
loading  level  approximately  horizontal. 

Test  chamber 

The  unconfined  compression  chamber 
in  the  hydraulic  press  (Fig.  9)  is  basically 
a  frame  with  a  leveling  base  upon  which 
the  test  specimen  rests.  The  loading  pis¬ 
ton  mounted  in  recirculating  ball  bushings 
provides  a  base  plate  for  the  load  measur¬ 
ing  transducer.  This  transducer  is  in 
direct  contact  with  the  top  spherical  sur¬ 
face  of  the  test  specimen  end  cap.  This 
arrangement  permits  measurement  of  the 
load  applied  to  the  test  specimen  at  any 
time.  Average  deformations  are  meas¬ 


ured  by  two  linear  motion  potentiometers  mounted  diametrically  opposite  each 
other  on  the  circumference  of  the  load  transducer. 


Load  and  deformation  measurements 

All  loads  applied  to  test  specimens  in  the  hydraulic  and  constant  stress 
presses  were  measured  with  Baldwin-Lima-Hamilton  load  cells  having  appro¬ 
priate  load  ranges. 

Hydraulic  press  loads  were  measured  using  the  load  cells  with  readout 
on  one  channel  of  a  Leeds  and  Northrop  Azar  G- type  X-X  recorder.  After 
calibration,  loads  were  measured  continuously  to  within  1.0%  of  the  applied 
load. 

Average  axial  deformations  of  test  specimens  in  the  hydraulic  press  were 
measured  using  two  carbon-strip,  infinite-resolution,  resistance-type  linear- 
motion  potentiometers  mounted  diametrically  opposite  each  other  on  the  load 
cell  (see  Fig.  9),  and  movements  were  recorded  on  one  channel  of  the  L  and  N 
recorder.  Using  calibration  charts,  deformations  were  measured  to  within 
0.  0025  in.  for  movements  less  than  0.25  in.  and  0.005  in.  for  movements 
greater  than  0.25  in. 
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Constant -stress -apparatus  loads  were  measured  with  the  load  cell  and 
lead  manually  using  a  B-L-H  type  N  portable  strain  indicator.  Loads  were 
determined  accurately  to  within  0.  3'fo  of  the  applied  load.  Deformations  were 
measured  using  dial  indicators  with  1/10,000  in.  gradations  and  a  sensitivity 
of  ii/J  00,  000  in.  (See  Figure  7  for  arrangement  of  load  cell  and  dial  indica¬ 
tors.  ) 

Constant  load  lever-type-press  loads  were  determined  by  computing  the 
hanger  weights  using  the  lever  arm  ratio.  The  load  applied  to  each  specimen 
was  checked  by  placing  a  load  cell  in  the  specimen  test  space  and  reading 
the  load  with  the  hanger  weights  in  place.  Deformations  were  measured 
using  the  same  type  of  extensometer  as  was  used  in  the  constant  stress  appa¬ 
ratus  (see  Fig.  7). 

Temperature  control 

Test  temperatures  of  25F  and  lower  in  the  walk-in  cold  room  were  con¬ 
trolled  to  within  ±!F.  To  damp  temperature  fluctuations  to  less  than  ±0.  hOF, 
tests  were  conducted  in  inclosures  constructed  of  Z  in.  thick  rigid  type  in¬ 
sulation  (Styrofoam). 

Test  temperatures  above  Z5F  were  controlled  by  heating  and  circulating 
air  within  the  insulated  test  inclosures.  Each  test  specimen  was  housed  in  a 
split  Lucite  cylinder  to  reduce  temperature  fluctuations  of  the  air  surrounding 
it.  Heat  was  supplied  by  light  bulbs  mounted  in  the  fan  air-stream.  The  tem¬ 
perature  was  regulated  by  a  mercury-column -type  thermoregulator  which 
activated  a  relay  to  supply  heat  upon  demand. 

Air  temperatures  within  the  Lucite  inclosures  surrounding  the  test  speci¬ 
men  were  held  constant  well  within  ±  0.  IF  of  the  desired  temperature. 

Temperature  measurements 

Temperatures  Z 5F  and  lower  were  measured  to  the  nearest  0.  sF  with  a 
mercury  thermometer  placed  within  the  insulated  inclosures.  Inside  the 
Lucite  inclosure  a  thermistor  sensed  the  temperature  of  the  air  surrounding 
the  test  specimen  and  the  readings  were  recorded  every  I.ZS  min  on  a  lii- 
point  L  and  N  type  H  recorder,  to  the  nearest  0.ZF.  Test  temperatures  above 
Z 5F  were  measured  to  within  0.  IF  using  the  same  sensing  and  recording  sy  - 
tern  but  with  increased  sensitivity.  Thermistor  readings  were  checked  daily 
using  a  manually  operated  Wheatstone  bridge.  Cold-room  temperatures  (out¬ 
side  of  .nclosures)  were  recorded  continuously. 


MATERIALS 

Ottawa  sand  (<i0-3t>)  and  Manchester  fine  sand  were  tested.  Ottawa  sand 
represents  a  nearly  idealized  granular  soil.  Manchester  fine  sand  i.,  a  natu¬ 
ral  sand  with  quite  uniform  gradation;  it  is  finer  than  Ottawa  sand  (see  Fig.  10 
for  gradation). 

A  limited  number  of  ice  specimens  frozen  under  the  same  conditions  as 
those  for  the  soils  were  tested  as  a  correlation  material  and  to  obtain  ice-test 
data  using  the  same  testing  equipment.  Ice-specimen  densities  were  less 
than  that  for  solid  ice  Decause  of  air  entrapment  during  freezing.  (See  Table 
AIII,  App.  A,  for  densities.  ) 
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Figure  1°.  Gradation  curve,  Manchester 
fine  sand. 


PREPARATION  AND  FREEZING  OF  TEST  SPECIMENS 
Specimen  molding 

The  soil  specimens  were  molded,  saturated  and  frozen  in  the  Lucite  mold 
previously  described.  To  reduce  friction  duri'  g  sample  ejection  from  the 
mold,  the  outsides  of  the  split  sleeves  were  coated  with  silicone  grease.  The 
soil  sample  was  not  in  direct  contact  with  the  grease.  However,  the  split  in 
the  forming  sleeve  allowed  saturating  water  to  contact  the  sil  cone  grease  thus 
permitting  only  a  slight  chance  of  silicone  contamination  of  the  specimen. 

After  the  mold  was  assembled  with  the  sleeves  in  position,  the  air-dry  sands 
were  placed  in  the  mold  in  six  equal  layers,  each  tamped  10  times  with  a  1/4 
in.  diam  steel  rod  (specimen  densities  are  tabulated  in  Tables  la,  Ila,  and 
Ilia).  With  top  and  bottom  of  the  mold  sealed,  the  soil-filled  mold  was 
evacuated  to  about  30  mm  of  mercury  using  a  water  ejection  pump,  and  then 
the  soil  was  saturated  from  the  bottom  by  admitting  de-aired  water  under 
vacuum  into  the  bottom  of  the  mold.  After  the  mold  was  filled  with  water 
an  additional  lOquarts  of  de-aired  water  percolated  through  the  specimens  to 
remove  as  much  trapped  air  as  possible.  When  saturation  was  completed 
the  top  and  bottom  mold  connections  were  sealed  for  placement  of  the  mold  in 
the  freezing  cabinet. 
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F  reezing 

After  saturation,  the  specimen-charged  moid  was  placed  in  the  freeing 
cabinet.  Spaces  between  the  sides  of  the  mold  and  cabinet  were  insulated  with 
granular  cork.  After  removal  of  the  top  mold  cover  a  de-aired  water  supply 
was  connected  to  the  bottom  of  the  mold  to  permit  specimen  freezing  in  an  open 
system.  In  this  arrangement  the  bottoms  of  the  specimens  were  exposed  to 
40F  with  a  free  water  supply,  and  the  tops  were  exposed  to  cold  circulating 
freezing  air.  The  rate  of  progress  of  the  32F  isotherm  was  determined  by 
means  of  thermocouples  spaced  1  in.  apart  along  the  vertical  axis  of  the  cen¬ 
ter  specimen.  Soil  samples  were  frozen  within  a  period  of  4  days,  i  e.  ,  the 
average  rate  of  frost  penetration  was  1  1/2  in.  per  day. 

In  freezing  ice  samples  the  procedure  was  similar  to  that  outlined  for 
soils.  However,  to  reduce  entrapment  of  air  bubbles  in  the  ice  the  rate  of 
freezing  was  reduced  to  approximately  1/4  in.  per  day. 

To  remove  the  frozen  specimens  from  the  freezing  mold,  the  mold  was 
clamped  in  a  specially  constructed  frame  and  each  specimen  was  ejected  by 
pressing  against  the  split  sleeve  using  a  hydraulic  jack  with  an  ejector  plate. 

The  ejector  plate  diameter  matched  the  mold  bore  and  had  a  raised  outer  rim 
that  permitted  ejection  of  the  specimen  in  the  sleeve  without  loading  the  speci¬ 
men.  The  split  sleeve  was  removed  from  the  ejected  specimen  by  expanding 
its  diameter  with  a  thin  blade  inserted  into  the  split. 

Trimming 

After  removal  of  the  split  sleeve,  the  specimen  was  inspected  for  imper¬ 
fections  and  cut  to  approximately  a  6-in.  height.  Rough  cutting  was  accom¬ 
plished  with  a  hacksaw  and  a  wooden  miter  box.  The  ends  were  squared  using  a 
special  case-hardened  vee-shaped  miter  box  and  various  gradations  of  wood 
rasps  and  steel  files.  After  the  specimen  ends  were  trimmed  to  final  length, 
the  dimensions  were  determined  by  measuring  the  circumference  (at  the  top, 
midheight,  and  bottom)  and  length  (at  six  locations  around  the  perimeter)  to  the 
nearest  1  /1000  in.  Variations  in  specimen  length  around  the  circumference 
were  within  ±  0.  005  in.  of  the  average  except  for  about  10  specimens  with  varia¬ 
tion  up  to  ±  0.  01  in.  of  the  average.  The  diameter  varied  less  than  ±  C.  003  in. 
along  the  specimen  length. 

The  volume  of  each  specimen  was  determined  by  submergence  in  liquid 
isooctane  (2,  2,  4  trimethylpentane)  at  20F. 

After  volume  determinations  were  made,  Ottawa  sand  sample  ends  were 
capped  with  a  thin  layer  of  ice  to  avoid  local  stress  concentrations  at  the  rela¬ 
tively  large  sand  grains  in  contact  with  the  loading  caps.  The  capping  was  ac¬ 
complished  by  pouring  a  layer  of  32F  water  on  a  flat  glass  plate,  then  setting 
the  specimen  end  on  the  plate  and  allowing  the  water  to  freeze  to  the  specimen 
in  the  cold  room.  The  glass  plate  was  removed  by  gently  warming  the  outer 
surface  of  the  plate. 

It  was  not  considered  necessary  to  cap  Manchester 
ice  since  the  fine  soil  grains  provided  a  smooth  contact 
caps. 

All  specimens  were  sealed  in  circumferential  rubber  membranes  and  metal 
end  caps  for  testing. 

Storage  and  tempeiing 

Prior  to  preparation  for  testing,  the  specimens  usually  remained  in  the 
sealed  freezing  mold.  Occasionally  it  was  necessary  to  eject  several  specimens 


fine  sand  specimens  with 
surface  for  the  loading 
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in  advance  of  preparation  for  testing.  These  specimens  were  sealed  in  rub¬ 
ber  membranes  and  temporary  etui  caps,  Lhen  sealed  in  plastic  bags  with 
crushed  ice  to  eliminate  sublimation.  No  change  in  specimen  weight  could 
be  detected  in  weighings  of  selected  specimens  before  and  after  storage. 
Storage  periods  did  not  exceed  6  weeks. 

Before  testing,  all  specimens  were  stored  at  the  test  temperature  for  a 
minimum  of  48  hours.  The  required  tempering  time  was  checked  using  three 
thermocouples  embedded  at  the  midpoint  and  quarter  points  of  the  axial  height 
of  a  control  specimen.  This  check  showed  tliaL  24  hours  was  sufficient  time 
for  the  sp'cimen  to  reach  equilibrium  at  the  test  temperature. 


CREEP  AND  STRENGTH  TESTING  PROCEDURE 

At  each  test  temperature  a  series  of  compression  type  tests  was  con¬ 
ducted  by  first  determining  the  "instantaneous"  strength*  of  the  frozen  soil 
or  ice  and  then  performing  creep  tests  at  reduced  stress  levels.  Each  test 
series  included  constant  stress  and  constant  load  compression  tests  per¬ 
formed  at  stress  levels  of  approximately  60,  35,  20,  10  and  2.5%  (stress 
ratio,  p)  of  the  average  instantaneous  strength.  One  Lest  series  was  conduc¬ 
ted  at  each  of  the  following  test  temperatures:  15,  25,  29  and  3  IF.  All  com¬ 
pression  type  tests  were  performed  on  specimens  2.8  in.  in  diameter  by  6  in. 
high.  Whether  constant-stress  or  constant-load  tests  were  performed  de¬ 
pended  upon  the  magnitude  of  the  applied  stress  and  the  expected  deformation. 
Constant  load  tests  were  used  for  high  stress  and  for  small  expected  defor¬ 
mation  at  low  stresses,  while  constant  stress  tests  were  performed  at  inter¬ 
mediate  stress  levels  (i.  e.  ,  in  the  range  of  p  -  15  to  40%)  where  the  deforma¬ 
tions  were  expected  to  be  large. 

The  compression  creep  test  on  each  specimen  was  performed  by  first 
applying  a  seating  load  of  approximately  2  psi  to  the  specimen  to  insure  posi¬ 
tive  contact  between  the  test  specimen  and  components  of  the  loading  system, 
and  then  applying  the  test  load  in  less  than  2  sec.  Instantaneous  strengths 
were  determined  by  moving  the  loading  press  head  against  the  specimen  at  a 
rate  of  0.2  in.  per  min  (strain  rate  of  0.  033  per  min).  After  each  test  was 
completed,  photographs  were  taken  of  the  lest  specimen  (see  Fig.  11-13  for 
typical  specimens)  and  wate  “  contents  were  determined. 


TEST  RESULTS 

Figures  14-16  show  typical  time-deformation  curves  (Ottawa  sand,  Man¬ 
chester  fine  sand  and  polycrystalline  ice)  produced  directly  from  data  for 
each  specimen  subjected  to  the  unconfined  compression  creep  test.  The  time- 
strain  curves  on  Figures  17  through  20  and  Figures  21  through  24  summarize 
the  data  for  Ottawa  sand  and  Manchester  fine  sand  respectively.  (Where  more 
than  one  specimen  number  is  shown  on  a  single,  curve,  the  curve  represents 
the  average  of  the  curves  obtained  for  specimens  indicated  and  the  vertical 
bars  indicate  the  total  range  of  values.)  A  summary  of  individual  specimen 
data  is  shown  in  Appendix  A. 


*  "Instantaneous"  strength  used  herein  is  the  maximum  strength  determined 
by  loading  the  test  specimen  at  a  constant  strain  ram  of  0.  03  3  per  min. 
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a.  Brittle  failure  in  instan¬ 
taneous  strength  test.  Temp, 
25F,  e~ 0.  58. 


Ows-51 

b.  Plastic-shear  failure  at 
815psi(63.5%of  instantaneous 
strength).  Temp,  25F,  e-0.60. 


c.  Specimen  subjected  to  200 
psi  (13.7%  instantaneous  strength) 
for  350o  hr  at  25F;  e  =  0,60. 


figure  it.  Typical  Ottawa  sand  specimens  after  testing. 
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Figure  2t.  Time  vs  strain,  Manchester  fine  sand,  31F. 


DISCUSSION 

Strain-time  data 

Creep  data  presented  as  deformation  vs  time  curves  are  the  primary  re¬ 
sults  of  this  investigation  and,  with  the  exception  of  the  "instantaneous"  stress- 
strain  curves,  are  the  basis  of  all  the  curves  reported  herein.  Typical  time- 
deformation  curves  for  sands  are  shown  in  Figures  14  and  15.  From  these 
curves,  true  ("logarithmic"  or  "natural")  strains*  were  computed  by  subtract¬ 
ing  test-machine  calibration  values  from  the  deformations  and  adjusting  the 
constant  load  test  values  to  a  constant  stress  basis.  Typical  time-strain  curves 
are  shown  in  Figures  17  -  24. 

Deformation  vs  time,  and  strain  vs  time,  curves  for  each  temperature  can 
be  grouped  according  to  shape.  The  curves  for  lower  stress  ratios  (<  15%)  how 
the  rate  'of  strain  decreasing  continuously  with  time  and  the  total  strain  (or  de¬ 
formation)  approaching  a  constant  value  asymptotically.  Intermediate  stress 
ratio  (15  to  60%)  curves  are  similar  in  shape  to  the  classical  creep  curves  shown 
in  Figure  1  and  exhibit  large  deformations  prior  to  failure,  f  Curves  for  high 
stress  ratios  (>  60%)  approach  a  straight  line  with  relatively  small  deforma¬ 
tions  prior  to  failure.  Figures  17  to  20  (Ottawa  sand)  and  21  to  24  (Manchester 
fine  sand)  are  sets  of  average  strain  vs  time  curves  showing  the  three  differ¬ 
ently  shaped  curves. 

The  low  stress  ratio  group  of  curves  consisting  of  instantaneous  elastic, 
time-dependent  elastic,  and  plastic  components  represents  strains  (or  defor¬ 
mations)  resulting  from  stresses  smaller  than  the  long-term  strength!  of  the 
frozen  soil.  The  various  deformation  components  are  labeled  on  the  deforma¬ 
tion  rebound  curve  in  Figure  25a.  The  total  maximum  long-term  strain  of  the 

*  True  strain  =  In  l/(l-ec)  where  the  conventional  strain  ec~&L/L0  =  axial 
deformation/original  length. 

t  See  page  1  for  definition  of  terms. 
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specimens  tested  did  not  exceed  0.05  and  0.  10  in.  /  in.  for  Ottawa  sand  and 
Manchester  fine  sand,  respectively.  Ten  rebound  tests  performed  on  frozen 
Ottawa  sand  indicate  that  the  total  elastic  comp  jnent  is  less  than  15%  of  the 
total  deformation  and  that  the  amount  of  instantaneous  elastic  deformation  is 
generally  less  than  time -dependent  elastic  deformations  (see  Table  I).  The 
percentage  ranges  of  total  deformations  shown  in  Table  I  indie  *e  that  ir¬ 
reversible  (plastic  type)  deformation  is  dominant  even  at  stress  levels  below 
the  long-term  strength. 
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Table  1.  Types  of  deformations  from  rebounded  creep  tests, 
frozen  Ottawa  and. 

Percent  of  Total  Deformatio  i 

Delayed  Plastic  or  Rebound 


T  emp 

(F) 

S'  ress 
(psi) 

Instant* 

rebound 

elastic 

rebound 

Total 

rebound 

plasto -viscous 
deform 

time  observed 

(h£) _ 

15 

1  70 

0  to  2. 5 

1 . 0  to  7 

1 .  3  to  9 

91  to  99 

1/2 

24.  5 

50 

0  .o  7 

6  to  10 

8  to  14 

86  to  92 

382 

29 

100 

0  to  2 

1  to  2 

1  to  4 

96  to  99 

1 

*  Deformation  observed  10  sec  after  removal  of  load. 

At  the  intermediate  stress  ratios,  the  strain  vs  time  curves  display  the 
characteristics  of  the  classical  creep  curves  for  metals.  These  curves  show 
instantaneous  elastic  and  delayed  elastic  (i.e.,  viscoelastic),  viscous,  visco¬ 
plastic  and  plastic  strains  and  eventual  failure.  In  Figure  25b  these  compo¬ 
nents  and  the  three  creep  stages  are  labeled  on  the  time  vs  strain  curve  for 
an  actual  test  specimen. 

Strain  components  and  strains  at  critical  points  are  summarized  in  Table 
II.  Section  A  of  summary  Table  II  shows:  (1)  that  all  instantaneous  strains 
recorded  (point  a,  Fig. 25b)  are  less  than  0.6%;  (2)  that  the  range  of  strains 
where  viscous  flow  begins  (point  v,  Fig. 25b)  is  0.  3  to  4.  7%  for  Ottawa  sand  and 
3  to  15%  for  Manchester  fine  sand;'  (3)  that  the  range  of  strains  at  initial  visco¬ 
plastic  flows  (point  f,  Fig. 25b,  start  of  tertiary  creep)  are  0.6  to  8.0  and  12  to 
28%  for  Ottawa  sand  and  Manchester  fine  sand,  respectively.  It  should  be  noted 
that  the  larger  strains  experienced  by  Manchester  fine  sand  under  equivalent 
conditions  could  be  explained  by  the  greater  void  ratio  of  these  specimens. 

The  magnitudes  of  each  of  the  strain  components  in  percent,  compared  in 
Section  B  of  Table  II,  do  not  reveal  an  obvious  trend  but  do  indicate  that  the 
vis  oelastic  component  (strain  component  3,  Figure  25b),  and  the  viscous 
component  (strain  component  4,  Figure  25b)  are  roughly  of  the  same  order 
of  magnitude  and  that  ihe  instantaneous  strain  is  insignificant  in  compari¬ 
son. 

Section  C  of  Table  II  compares  the  three  strain  components  as  percentages 
of  the  total  strain  experienced  by  test  specimens  prior  to  the  nset  of  final  visco¬ 
plastic  strain  (i.e.,  as  percent  of  the  strain  at  point  f).  At  each  temperature  in 
this  tabulation,  the  viscous  component  of  the  Manchester  fine  sand  curve  tends 
to  increase  and,  correspondingly,  the  initial  viscoelastic  components  tend  to 
decrease  as  the  stress  increases.  This  trend  does  not  necessarily  mean  that 
the  viscoelastic  component  is  smaller  but  that  it  can  be  obscured  by  the  more 
stress-sensitive  viscous  component.  Rebound  creep  tests  at  the  higher  stress 
levels  would  help  clarify  this  point.  Plots  showing  these  trends  are  given  in 
Figures  26a  and  26b.  No  similar  trends  were  detected  in  the  Ottawa  sand  data, 
possibly  because  the  smaller  void  ratio  of  the  Ottawa  sand  permits  less  viscous 
flow  of  the  ice  matrix  in  the  sand  voids  prior  to  the  development  of  greater 
frictional  resistance  between  the  sand  grains. 

It  was  noted  for  stresses  greater  than  the  long-term  strength  that  viscous 
type  deformation  was  observed  earlv  in  the  tests,  an  observation  that  permits 
an  early  prediction  of  eventual  specimen  incipient  failure  and  failure.  *  Failure 
occurred  by  either  an  abrupt  brittle-type  fracture  (Fig.  11a)  or  a  plastic  flow 
accompanied  by  fissuring  (Fig.  lib  and  12b).  The  plastic-fiow  type  failure 
occurred  after  large  deformations,  and  fissures  often  preceded  the  abrupt 
collapse  of  the  specimen. 


*  See  p.  I  for  definition  of  terms. 


Table  II.  Strain  component*  and  strain  at  critical  points  under  intermediate  stresses  (see  Fig.  2  5b). 
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High  stress  ratio  curves  exhibit  an  abbreviated  form  of  the  classical  creep 
curve  passing  rapidly  through  the  three  stages  of  the  creep  curve  and  terminat¬ 
ing  in  an  abrupt  failure. 

It  should  be  noted  that  the  deformation  curves  for  ice  (Fig.  16a,  b)  exhibit 
the  same  general  types  of  curves  corresponding  to  each  of  the  three  stress 
ratios  for  the  frozen  sands  except  that  low  stress  ratios  in  ice  produced  un¬ 
damped  viscous  deformation  (i.  e.  ,  a  constant  rate  of  deformation)  rather  than 
damped  deformation  exhib' ted  by  the  sand.  Clen  (1955),  Jellitiek  anu  Brill 
(1956)  and  Butkovic'n  and  Landauer  (1959,  1960)  investigated  creep  and  visco¬ 
elastic  properties  of  ice  under  low  stress  and  have  formulated  theories  and 
equations  describing  its  behavior  as  influenced  oy  stress,  temperature  and 
structure.  Tests  on  ice  in  the  present  investigation  were  conducted  solely  for 
comparison  with  similar  tests  on  frozen  sand  under  the  same  test  conditions. 

Strain  rate 

P.ates  of  strain  were  determined  as  the  slopes  of  the  strain  vs  time  curves 
using  a  digital  computer.*  The  rates  for  each  test  specimen  were  computed  by 
considering  successive  groups  of  five  consecutive  data  points  on  the  strain  vs 
time  curves.  A  second  degree  polynomial  equation  was  fitted  to  the  five  data 
points  by  the  method  of  least  squares,  then  the  slope  of  the  equation  was  deter¬ 
mined  at  the  middle  point  of  the  five  data  points.  New  groups  of  five  data  points 
were  considered  by  advancing  along  the  strain  vs  time  curve  one  point  at  a 
time  (i.e.,  by  eliminating  the  first  data  point  and  including  a  new  advanced  data 
point).  The  process  of  fitting  the  polynomial  equation  and  determining  the  slope 
of  the  equation  was  repeated  for  each  group  of  five  data  points.  Using  this  sys¬ 
tem,  rates  of  strain  were  determined  for  the  entire  length  of  the  strain  vs  time 
curve  except  for  the  first  tv  o  and  last  two  points  on  the  curve.  The  procedure 
was  readily  adapted  to  digital  computer  methods  and  by  considering  only  five 
points  at  a  time  it  was  not  necessary  to  determine  an  expression  for  the  entire 
strain-time  curve. 

Typical  strain  rate  vs  time  curves  for  both  the  damped  and  undamped  type 
creep  tests  are  plotted  on  Figure  27.  The  undamped  curve  has  a  high  initial 
strain  rate  which  decreases  hyperbolically  (during  creep  stage  I)  to  a  minimum 
relatively  constant  value  (creep  stage  II  or  viscous  flow,),  then  increases  rapidly 
(creep  stage  111)  until  specimen  failure.  The  initial  portion  of  the  damped  curve 
is  similar  in  shape  to  the  undamped  curve  except  that  it  approaches  a  zero 
strain  rate. 

Typical  3train-rate  vs  reciprocal-of-time  curves  for  three  applied  stresses 
plotted  on  log  coordinates  (Fig.  28)  show  that  the  strain  rates  of  both  the  damped 
and  undamped  creep  curves  are  initially  straight  lines.  The  undamped  curve 
represented  by  the  260  psi  curve  reaches  a  minimum  strain  rate  then  increases 
rapidly  as  indicated  by  the  sharp  turn  to  the  right.  In  contrast  the  damped  curve 
remains  a  straight  line  then  decreases  (i.  e.  ,  turns  to  the  left)  as  the  rate  of 
strain  approaches  zero  (see  Fig.  29).  Although  the  20  psi  curve  on  Figure  28 
is  for  a  damped  creep  test,  it  does  not  show  the  decrease  in  creep  rate  as 
clearly  as  the  damped  curves  on  Figure  29.  One  possible  explanation  for 
this  decrease  in  strain  rate  can  be  that  friction  between  sand  grains  increases 
enough  to  dominate  the  creep  of  the  test  specimen  and  eventually  produce 
stability. 


*  The  method  of  determining  strain  rates  was  suggested  by  Dr.  A.  Assur, 
Chief  Scientist,  USA  CRREL. 
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'igure  28.  Creep  rate  and  reciprocal  of  time,  Ottawa 
sand,  3 IF. 


Figure  29.  Creep  rate  and  reciprocal  of  time,  Man¬ 
chester  fine  sand,  15F. 


The  three  damped  curves  on  Figure  29  for  the  constant  temperature  of  1  5F 
indicate  that  not  only  the  magnitude  of  the  strain  rate  but  also  rate  of  change  of 
the  strain  rate  (i.e.,  the  slope  of  the  curve)-is  stress  dependent.  Since  the 
slope  is  steeper  for  higher  stresses,  the  creep  decelerates  faster  and  stabilizes 
faster  for  the  higher  damping  stresses. 

By  comparing  the  )  5F  and  25F  curves  of  Figure  30,  it  is  clear  that  at  in¬ 
creased  temperatures  the  rate  of  strain  is  increased.  The  10  psi  difference  in 
test  stresses  between  the  two  curves  produces  a  relatively  minor  cltange  at  the 
temperatures  involved.  A  similar  comparison  between  the  29F  and  31F  curves 
(Fig.  30)  shows  that  the  rate  of  decrease  in  strain  rate  is  greater  at  the  higher 
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Figure  30.  Creep  rate  and  reciprocal  of  time,  Manches¬ 
ter  fine  sand  at  various  temperatures. 


test  temperature  even  though  the  stress  is  smaller  (i.e.,  the  slope  is  greater  at 
the  higher  temperature).  This  indication  of  deceleration  of  creep  with  increased 
temperature  requires  further  investigation  since  only  a  trend  was  noticed  here. 


Strain  equations 

An  equation  to  be  used  for  the  prediction  of  the  creep  strain  of  a  given  frozen 
soil  must  at  least  include  factors  that  reflect  the  influence  of  the  soil  tempera¬ 
ture,  the  magnitude  of  the  applied  stress,  the  time  during  which  the  stress  acts, 
the  ice  or  water  content,  the  soil  particle  size  and  the  density.  Although  other 
factors  affect  the  creep  strain,  their  influence  is  reflected  indirectly  by  these 
factors  listed. 

Strain  equation  by  Vialov.  S.  S.  Vialov  (1962)  has  presented  a  creep  strain 
equation  based  on  tests  performed  on  undisturbed  sandy  silt  and  clay,  which  in¬ 
cludes  the  most  essential  factors  either  directly  or  indirectly.  Vialov's  equa¬ 
tion  for  total  strain  is: 

r  o-t*  I  l/m 

£  (0  +  l)k  J  +■  *o  (!) 


where:  cr 

t 
6 

X,  m,  u  and  & 


applied  constant  stress  kg/cmz 
time,  hr 

temperature  in  degrees  C  below  the  freezing  point  of  water 
instantaneous  strain 

constants  that  are  characteristic  of  the  material  (and 
depends  on  units). 


This  equation  conceals  the  fact  that  in(l  +  9),  the  1  has  temperature  units.  For 
this  reason,  Assur  (1963)  suggested  that  the  equation  be  written: 
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Here,  the  numerical  values  of  m  and  X  do  not,  but  ui  and  90  do,  vary  with  the 
units  employed.  90  =  an  assumed  constant  reference  temperature  greater 
than  zero.  Vialov  assumed  9e  =  1C  so  that  lg(H>0)  had  meaning  at  0  =  0.  In 
this  report  90  =  IF  was  used.  For  numerical  comparisons,  stress,  time  and 
temperature  data  should  be  converted  to  the  same  units. 

The  first  term  on  the  right  side  of  this  equation  was  developed  from  the 
power  function  relationship,  a  =  A  tm,  b  .tween  stress  and  strain  for  stresses 
less  that  the  long-term  strength  of  the  test  material,  where  A  =  f(9,  t),  in  the 
units  of  stress. 

Values  for  A  and  m  are  obtained  from  log  plots  of  stress  against  strain 
(see  Fig.  31a.  b)  for  various  time  intervals  after  stress  application  to  the 
test  specimens.  Vialov  represented  sandy  silt  and  clay  data  by  a  straight  line 
on  a  log  (rvs  log  «  plot  even  though  he  recognized  that  a  curve  did  exist.  The 
sand  data  reported  herein  also  showed  some  curvature.  This  curvature  along 
with  the  experimental  scatter  made  it  difficult  to  represent  the  log  strain  vs 
log  stress  curves  as  straight  lines.  However,  to  determine  the  constants  in 
Vialov's  strain  equation  for  the  sand  data,  straight  lines  were  drawn  as  indi¬ 
cated  in  Figures  31a  and  b  for  all  test  data.  The  average  slopes  of  these 
curves  (m)  are  shown  in  Table  III.  This  table  indicates  that  within  the  accuracy 
of  the  tests,  the  values  of  m  are  independent  of  temperature  for  both  the  Ottawa 
sand  and  the  Manchester  fine  sand,  m  may  be  considered  a  characteristic  of 
the  material  and  its  density,  and  its  value  is  independent  of  the  units  employed. 


Tablo  HI.  Value  of  m  in  Vialov's  strain  equation. 


Temperature 
°  F 

Ottawa 

sand 

Manchester 
•fine  sand 

15 

0.  78 

0.  42 

25 

0.79 

0.  36 

29 

0.  77 

0.  40 

31 

0.79 

0.  35 

Average 

0.78 

0.  38 

Table  IV. 


Constants  for  Vialov's  strain  equation. 

u  k 


Material 


For  90  =1 
si(hr)X]  /*Fk 


For  90  =  1*C 
si(hr)Xl  /*Fk 


Ottawa  sand  (20-30  mesh) 

0.  78 

0.  35 

5500 

0.  97 

3600(456)t 

1.0 

Manchester  fin:  sand 

0.  38 

0.  24 

285 

0.  97 

185(23.4) 

1.0 

(40-200  mesh) 

Silt*  (Callovian  sandy  loam) 

0.  27 

0.  10 

90 

0.  89 

76(  9.0) 

0.  89 

Clay*  (Bat-Baioss  clay) 

0.  4 

0.  18 

130 

0.  97 

103(12.8) 

0.  97 

*  Data  from  Vialov  et  al_.  (1962),  Chapter  V. 
t  oj  in  . kg/cm2  hrk/*Ck  shown  in  parentheses. 
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By  plotting  A  (i.e.,  the  ordinates 
for  log  strain  equal  to  unity)  from  Fig¬ 
ures  31a  and  b  against  time  for  the  dif¬ 
ferent  test  temperatures,  values  for 
factor  R(=  At^-)  and  exponent  X  were 
obtained  as  shown  on  Figure  32.  The 
constants  to  and  k  were  similarly  ob¬ 
tained  by  plotting  the  intercept  R 

Q 

against  9o(tt~  *'  0  as  shown  on  Figure 

3  3.  The  values  for  Ottawa  sand  and 
Manchester  fine  sand  are  listed  in 
Table  IV  along  with  those  for  sandy 
silt  and  clay  published  by  Vialov  (1962). 


/ /  -  Figures  34a  and  b  show  a  compar- 

/.////  _  ison  between  Vialov's  equation  and  ac- 

tual  test  data  for  Manchester  fine  sand 
10  —  and  Ottawa  sand,  respectively. 

Strain  equation  based  on  strain  rate. 
Vialov's  equation  demands  quite  pre- 
I  i  I  1 1 . 1.!  I  ■  1  ■  1 1  hi  I  i  I  1 1  .i.l  cise  measurements  of  the  absolute  val- 
1  10  l0*  ues  of  strain  (i.e.,  deformation)  and 

temperaiure.  e  also  requires  several  tests  for  the 

evaluation  of  the  formula  constants. 

Figure  33.  Factor  R  and  temper  ature.  To  overcome  these  two  difficulties  and 

to  provide  a  more  flexible  approach,  a 
strain  equation  was  developed  using 

strain  rates  where  only  the  differences  in  strain  (i.e.,  deformation)  need  to  be 
measured.  Like  Vialov's  equation  the  only  stress  condition  under  considera¬ 
tion  is  where  the  applied  stress  is  less  than  the  long-term  strength  (i.e.,  for 
damped  creep  tests)  and  the  instantaneous  strain  is  neglected. 

The  straight  line  portion  of  the  log  (1/t)  vs  log  strain  rate  curves  for  damped 
creep  tests  on  Figures  28  to  30  can  be  represented  by: 


Figure  33.  Factor  R  and  temper  ature. 


-1/M 


By  integration 


■  *.(*). 


(4r) 


where:  €  =  strain 

«  =  strain  rate 

ij  =  strain  rate  1  hour  after  stress  is  applied 
e  o  =  initial  instantaneous  strain  as  stress  is  applied 
t  =  time  after  stress  is  applied 
M  -  slope  of  log  1/t  vs  log  strain  rate. 

An  examination  of  the  strain  rate  curves  (Fig.  30)  reveals  that  M  is  a  func¬ 
tion  of  both  temperature  and  stress,  but  for  the  temperature  range  investigated, 
the  stress  effect  predominates.  Assuming  M  is  independent  of  temperature,  and 
is  a  straight  line  function  of  stress  on  logarithmic  coordinates  (Fig.  35),  the 
equation  for  M  becomes: 

M  -  C  <r1/w  (3) 


r 
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b  Ottawa  SAND  {20-30),  170  pit 

Figure  36.  Strain  and  time,  15F. 

* 

where  w  is  the  slope  of  log  o -  vs  log  M  curve  and  treated  as  a  constant  for  each 
material.  Values  of  w  determined  graphically  for  the  Ottawa  sand  and  Manches¬ 
ter  fine  sand  at  the  densities  tested  are  9  and  35,  respectively,  and  C  is  almost 
exactly  unity  for  each  sand. 

The  values  cf  «  .  and  M  in  eq  2  can  be  determined  graphically  using  the 
slopes  of  the  cangent  to  the  strain  vs  time  curve  on  arithmetic  coordinates  at 
times  of  |  hour  and  1  hour  after  stress  application  (Fig.  36a,  b)  or  by  using 
the  log  t  vs  log  1  it  curves  for  the  first  8  hours  of  the  creep  test  data  (see 
Fig.  29).  The  arithmetic  coordinate  method  offers  a  simple,  rapid  means  for 
predicting  strain  but  is  not  as  accurate  as  the  method  using  the  log  «  vs  log  1  / 1 
curves.  The  logarithmic  coordinate  method  predicts  long-term  strains  that 
are  in  close  agreement  with  the  actual  data;  however,  the  method  requires  the 
determination  of  the  strain  rates  at  several  points  on  the  8-hour  time-strain 
curve.  These  strain  rates  can  be  found  graphically  by  drawing  tangents  as 
indicated  in  Figures  36a  and  b  or  by  fitting  a-polynomial  expression  to  a  portion 
of  the  strain-time  curve  and  then'differentiating  it,  a  method  that  was  used  to 
obtain  the  strain  rates  in  Figures  29  and  30.  Strain  curves  predicted  by  eq  2 
using  both  arithmetic  coordinates  and  logarithmic  plotting  for  determining  «j 
and  M  are  compared  with  test  data  on  Figures  34a  and  b. 

A  more  general  expression  for  can  be  obtained  by  representing  the 
logarithmic  plots  of  tj  vs  <r  for  each  temperature  by  a  straight  line  as  shown 
in  Figures  37a  and  b.  The  general  equation  for  each  temperature  is: 


r.  STRESS. 
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uf»o*a  .x 
SAND/ 


€,.  STRAIN  RATE  AT  ONE  HOUR,  hr 
b  Manchttfor  fin*  SAND 

Figure  37.  Strain  rate  at  time  1  hr,  and 
stress.  Various  temperatures. 


+  0  temperature  ,'f 

Figur  38.  Temperature  and  stres s  o-i 
for  unit  strain  rate  at  I  hr  ( Vj ) . 


=  ffl  *1 


i,  =  (f-)' 

1  <*i 


where 


<Tj  =  stress  at  which  J  is  unity 
K  =  constant  slope  of  the  straight  line  plot. 


This  equation  assumes  that  o"i  is  temperature  dependent,  and  that  K  is  con¬ 
stant  for  the  materials  and  the  temperature  range  investigated  (as  indicated  by 
the  parallelism  of  the  lines). 

Temperature  and  <r;  are  plotted  on  Figure  38.  [The  quantity  { W  +-1 )  =  Q0  (£—  +  1) 

i  0 

was  used  instead  of  Q  to  avoid  the  logarithm  of  zero  at  a  temperature  of  32F.  ] 

The  equation  for  each  straight  line  on  Figure  38  is: 

<fi  =  <%i  (Q  +  1)Q  (5) 

where  (%i  is  stress  at  9  =  0  (temp  32F) 

q  (const.)  is  slope  of  the  straight  line  plot. 

By  substituting  the  expression  for  <Tj  into  eq  4,  an  expression  for  J  is: 


(6) 
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‘  [ff0i  (0+1)  J 

where  9  is  in  Fahrenheit  degrees. 

Using  '■  eq  2  becomes 


*i  Oj  +  *o 


and  substituting  j  from  eq  6,  then 

1  /  K  4< 

_ 1 _  + 

nQ/  ®  ,  nQ  4<  * 0 

<5i  9  (y  +  1) 


o-qi  (9  +  1) 


t4- 

—  + 


On  Figures  34a  and  b,  eq  2a  and  b  (neglecting  (())  are  compared  with  Via¬ 
lov's  equation  and  the  test  data.  Equation  2a  has  the  advantage  of  requiring 
only  a  single  short-term  creep  test  to  predict  the  total  strain  for  a  given  mat¬ 
erial  after  applying  a  stress  for  a  specified  time.  Equation  2b  can  be  used  to 
estimate  total  strain  where  only  the  general  type  soil  (i.e.,  coarse  sand,  fine 
sand,  silt,  clay,  etc.)  is  known  and  constants  for  the  soil  type  can  be  estimated 
by  comparing  the  soil  under  consideration  with  similar  soils  that  have  been 
tested  previously.  (Table  Va  lists  the  values  of  constants  for  the  sands.)  Of 
course,  if  time  and  facilities  are  available,  creep  tests  should  be  made  on  the 
actual  soil  under  consideration  to  determine  both  the  long-term  strength  and 
the  total  strain  for  any  particular  stress. 

Table  Va.  Constants  for  strain  equation. 


_ Material _ 

Ottawa  sand 
Manchester  fine  sand 


(9  in  F  deg) 


1  5,  000 
16,  400 


Table  Vb.  Comparison  of  constants  for  strain  equations. 


Material 

a 

K 

k 

m 

a/K 

k/m 

Ottawa  sand 
Manchester  fine  sand 

0.  58 

0.  46 

0.  76 
0.71 

0.  97 

0.  97 

0.  78 

0.  38 

0.  76 

0.  65 

1 . 24 
2.  55 

A  temperature  dependency  comparison  can  be  made  between  Vialov's 
modified  strain  equation, 

\  1/m 

o'  ^ 

v  -  — TT3-— T  +  < 

weoV5”+  1) 

and  the  strain  rate  equation, 


„q,9  , .a 

<%i  90  (q^-  +  1) 


t*  + 


By  neglecting  (  and  considering  strain  to  be  a  function  of  temperature  only, 
the  strain  functions  become: 
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Me)  =  fe0<Tj-  +  i)lk/n' 

M0)  =  fe0(§-  +  D1q/k 

y0 

and  for  90  =  1°  in  the  same  system  of  units; 
fv(0)  =  (9  +  l)k/m 

fs(0)  =  (e  +  i)a/K 

which  indicate  that  the  exponent  a/K  should  be  equal  to  k/m.  By  using  values 
given  in  Tables  IV  and  Va  and  by  making  the  comparison  shown  in  Table  Vb, 
it  is  clear  that  the  exponents  are  not  equal.  However,  if  we  substitute  test 
values  in»o  both  equations,  we  find  that  for:  Ottawa  sand  at  15F  (0  =  17F)  and 
it  ~  170  psi  eq  1  yields 


170  t0>  35 


5, 500(17  +  1)°- 97 


=  3.  3  x  10*4  t0-4' 


1/0. 78 


and  for:  w  =  9,  M  =  oM  w  _  170^/9  -1,77 
0/  = 


M  *  1  -  -L-7-1_eJ-  =  0.44 


TT 


T77T 


in  pii,  hr,  and  F  deg. 


equation  Zb  gives 

-  % 

*  fo 


170 


n  1/0. 76 


1  5,  000  (17  +  l}°-  58 
=  6.  85  x  10~4  t°-  44 


t0.  44 

-q  ■  in  psi,  hr  and  F  deg. 


which  shows  that  equation  1  gives  values  of  about  one  half  of  equation  Z.  Per¬ 
haps  M  is  not  entirely  independent  of  temperature  ;  t hi s  should  be  investigated 
further. 

Strength  -  time 

Using  the  time  to  total  failure  and  the  corresponding  applied  failure  stress 
from  the  strain  vs  time  curves,  curves  for  strength  vs  time  at  different  test  tem¬ 
peratures  were  plotted  on  Figures  39  and  40.  The  curves  are  drawn  to  approach 
asymptotically  the  value  of  the  maximum  applied  test  stresses  that  did  not  cause 
failure  during  the  test  period  (i.e.,  a  stress  known  to  be  less  than  the  long-term 
strength).  At  the  termination  of  each  test  in  which  the  specimen  did  not  fail, 
the  measured  rate  of  strain  was  either  zero  (i.e.,  no  deformation  change  could 
be  detected  for  at  least  100  hours),  or  the  rate  of  deformation  was  decreasing 
continuously  after  a  total  test  period  of  over  Z000  hours.  The  curves  and  test 
data  in  Table  VI  show  that  the  saturated  frozen  sands  have  long-term  strengths 
of  less  than  15%  of  their  instantaneous  strength  it  should  be  noted  on  Figures 
39  and  40  that  ice  strengths  are  less  than  those  of  the  sands  tested  under  the 
same  conditions. 

Vialov  (1959)  suggests  that  the  variation  of  strength  of  frozen  soils  with 
time  is  an  exponential  function  in  the  form 

<rult  =  ln(t/B') 

'here  (5  and  B  are  parameters  which  depend  upon  the  soil  and  its  temperature 
and  t  and  <ruit  are  the  time  to  either  failure  or  incipient  failure,  and  failure 
stress,  respectively.  For  convenience  in  graphic  plotting,  the  logarithm  to  the 
base  ten  is  used  to  evaluate  (3  and  B;  then  eq  7  becomes: 
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Table  VI.  I.ong-term  unconfined  compressive  strength. 


Temp 

(°  F ) 

Instantaneous 
strength 
(P  si) 

Ottawa  Sand  (20-30) 

Predicted 

Ion  g  -term 
strength 
(psi) 

C  reep  test 

Non-failure 

(psi) 

stresse 

Failure 

(PsiI_ 

i*  nf^ 

1  5 

2530 

420  572 

1  70 

600 

24 . 5 

1460 

254  263 

200 

460 

29 

1  320 

134  149 

1  00 

250 

31 

745 

69  76 

40 

80 

Manchester  Fine  Sand 

1  5 

2790 

329  343 

350 

560 

25 

2000 

176  193 

1  60 

400 

29 

1340 

115  128 

1  00 

265 

31 

805 

63  66 

80 

1  50 

*  I  -  predicted  long-term  strengths  are  for  100  years  determined  by  eq  7ausing 
Figures  41  ana  42  to  determine  values  for  B  and  0.  See  Table  VII. 

|II  -  predicted  long-term  strengths  are  the  100-year  strength  determined  by 
eq  7a  using  two  short -term  creep  strength  values  (e.p.  ,  60  and  35%  of  instan¬ 
taneous  strength)  to  evaluate  B  and  0. 

fable  VII.  Constants  foreq7aas  determined  from  Figures  41  and  42. 


O'tawa  Sand  Manchester  Fine  Sand 


Temp 

(°F) 

B 

(  hi') 

0 

(psi) 

B 

Uni _ 

0 

_ (Psi) _ 

1  5 

0. 0031 

3590 

0.  0126 

2600 

25 

0.  01 89 

1960 

0.  042  ! 

1290 

29 

0.  1 

1  000 

0. 0394 

852 

31 

0. 0994 

475 

193 

420 

0 

(7  a) 

^ult 

Togvt/BT 

The  values  of  parameters  0  and  B  can  be  determined  by  either  developing 
a  plot  of  1  /tr(  , .  vs  time  or  by  using  results  from  short-term  creep  tests. 

Figure's  4i  and  4 2  indicate  Vialov's  method  of  determining  these  parameters 
from  the  1/ir  ^  vs  log  time  curves.  These  curves  are  based  on  results  of 
both  short-term  and  long-term  creep  tests  and  can  be  used  to  predict  failure 
strengths  it  ore  accurately  than  by  using  only  short-term  creep  test  results 
alone.  The  disadvantage  of  using  these  plots  is  the  length  of  time  and  care 
required  to  perform  the  long-term  creep  tests.  Using  the  results  of  short¬ 
term  creep  tests  for  predicting  long  term  strength  has  the  advantage  of  re¬ 
quiring  only  two  simple  creep  tests  that  car  be  {  ermrmed  at  two  different  stress 
levels  within  a  period  of  3  days.  Values  of  0  and  B  can  be  determined  by  solving 
eq  7ausingthe  results  ol  these  two  tests.  For  comparison,  values  of  the  100- 
year  predicte  1  strengths  are  shown  ir.  Table  VI  in  columns  J  and  II  as  deter¬ 
mined  bv  the  >  / a u 1 1  vs  io^  time  curves  and  by  two  short-term  creep  tests,  res¬ 
pectively  . 
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ll  should  bo  noted  that  Vialov's  strength  equation  (eq  7)  should  not  be 
used  tor  extremely  short  periods,  since  when  the  duration  of  loading  ap¬ 
proaches  zero  it  predicts  strengths  approaching  infinity.  This  restriction 
does  not  invalidate  the  equation's  usefulness  in  predicting  the  long-term 
strength  of  the  frozen  soil. 

Further  examination  of  Table  VI  shows  that,  with  the  exception  of 
Manchester  fine  sand  at  i5  and  3  1  F ,  the  predicted  long-term  st  rength  for  both 
soils  at  the  different  temperatures  is  g  re  ate  r  than  the  tes  ted  long  -  te  rm 
strength  (non-failure  value  in  the  table)  and  less  than  the  next  higher  test  stress 
level  that  caused  failure.  Strength  variations  with  time  predicted  b^  eq  7a 
using  two  short-term  creep  tests  are  compared  with  test  strength  values  in 
Figures  30  and  40.  The  figures  and  Table  VI  show  that  the  Vialov  equation 
may  predict  unsafe  strength  values  if  based  on  short-term  creep  par¬ 
ameters  unless  an  appropriate  factor  of  safety  is  used. 

An  examination  of  Figures  39  and  4U  eveals  that  tl  e  greatest  loss  of 
strength  for  frozen  sand  occurs  during  a  relatively  shoit  period  after  applica¬ 
tion  of  stress.  By  comparing  values  in  Table  VIII,  it  is  clear  that  the  ability 
of  the  frozen  sand  to  resist  complete  failure  is  reduced  by  at  least  50'70  of  its 
instantaneous  strength  when  the  applied  stress  acts  24  hours  and  that  further 
strength  reduc'ions  occur  at  a  much  lower  rate.  Also,  stresses  that  are  to  be 
resisted  for  1000  hours  must  be  reduced  to  less  than  25'/o  of  the  instantaneous 
strength  of  the  frozen  sand.  It  should  be  emphasized  that  percentages  shown  in 
Table  VIII  are  to  demonstrate  the  strength  reduction  of  frozen  sand  and  are 
based  on  an  "instantaneous"  strength  defined  in  this  report. 


Table  VIII.  Percent  of  instantaneous  strength  loss  after  application  of  stress. 

Ottawa  Sand  (20-30) 

Time  after  stress  application 

Temp 


0  F 

1  hr 

24  hr 

100  hr 

1 000  hr 

1  5 

45 

55 

60 

75 

2  5 

25 

50 

55 

80 

29 

25 

50 

70 

80 

31 

45 

70 

80 

90 

Manchester  Fine  Sard 

1  5 

45 

65 

70 

80 

25 

4  5 

70 

80 

90 

29 

55 

75 

80 

95 

3  1 

50 

70 

80 

90 

Strength  -  temperature 

Except  for  the  instantaneous  strength  curve,  the  strength  vs  temperature 
curves  on  Figures  43  and  44  are  replottcd  from  the  strength  vs  time  curves 
for  the  various  times  of  stress  application  shown.  The  instantaneous  strength 
curves  are  drawn  through  the  average  instantaneous  strength  points  and  for 
comparison  purposes  instantaneous  ice  strengths  are  plotted  on  the  figures. 
The  long-term  strength  curves  arc  based  on  extrapolations  of  the  strength  vs 
time  curves. 
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TEMPERATURE  ,  *F 

Figure  43.  Strength  fur  various  conditions,  Ot¬ 
tawa  sand  (20-30). 


Figure  44.  Strength  for  various  conditions,  Man¬ 
chester  fine  sand. 
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In  Figure  4  3,  the  instantaneous  strength  curve  for  Ottawa  sand  shows  unex¬ 
pected  change  in  strength  values  near  29F.  A  similar  change  in  strength  in  this 
temperature  range  has  been  observed  by  others  (ACFFL,  1952).  It  was  noted 
that  at  25F  and  15F  brittle  type's  of  failures  occurred  and  at  31F  all  failures 
were  plastic.  One  possible  explanation  for  strength  irregularity  is  that  a  tran¬ 
sition  between  brittle  and  plastic  type  failure  occurs  at  about  29F.  The  type  of 
specimen  failure  for  a  constant  temperature  depends  upon  both,  density  and  rate 
of  loading.  The  influence  of  specimen  density  is  suggested  by  comparing  strength 
vs  temperature  curves  for  Manchester  fine  sand  reported  by  ACFEL  (1952)  with 
the  instantaneous  strength  curve  on  Figure.  44.  The  ACFFL  tests  performed  on 
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Figure  45.  Temperature,  |i  and  B  parameters. 

@0  =  IF 

specimens  with  void  ratios  of  0.  658  show  an  irregular  strength  curve  near  25F 
while  the  tests  shown  on  Figure  44  (average  void  ratios  of  0.  78)  did  not  reveal 
this  irregularity  even  though  the  rate  of  loading  was  over  five  times  greater 
than  that  applied  by  ACFEL.  Also,  all  specimens  reported  on  Figure  44  failed 
plastically.  A  further  investigation  of  instantaneous  strength  is  required  for  a 
positive  explanation  of  the  results  shown  in  this  report. 

The  curves  in  Figures  43  and  44  show  that  frozen  sand  strength  increases 
with  decreasing  temperatures  and  that  the  increase  in  long-term  strength  is 
less  for  a  given  temperature  reduction  than  the  short-term  strength  increase. 
For  example,  Ottawa  sand  showed  an  instantaneous  strength  increase  of 
<1000  psi  while  the  long-term  strength  increased  only  150  psi  for  a  tempera¬ 
ture  decrease  from  31F  to  15F. 
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For  the  temperature  range  tested,  the  long-term  strength  increase  with 
decreasing  temperature  may  be  approximated  by  the  equation  (rujt  -  a  +  b  |  0 ] 
(suggested  by  Vialov,  1959)  where  n,  a  and  b  are  constants,  b  has  the  units  of 
stress /degrees11. 

The  long-term  strength  of  Ottawa  sand  may  be  represented  by: 

a  h  =  20  ]  flj  in  psi  and  F  deg 
and  similarly  that  for  Manchester  fine  sand  may  be  expressed  by: 

3/4 

oujj.  =  40  +  40  |  0 ]  in  psi  and  F  deg. 

Figures  41  and  42  and  Table  VII  show  that  p  and  B  are  functions  of  tem¬ 
perature.  Plots  in  logarithmic  coordinates  of  p  and  B  vs  a  temperature  func¬ 
tion  are  shown  in  Figures  45a  and  b.  Equations  for  the  straight  lines  approxi¬ 
mating  the  plotted  data  for  Ottawa  sand  arr : 

0  q  1  Q  0.91 

p  =  26O0o°-91(1  +J-) 

Ho 

- 1  6  9  •”  *  ’  ^ 

B  =  0.  64  90  (i  +“) 

ao 

For  90  =  IF 

p  =  260(1  +  0)°-91 
B  =  0.64(1  +  0)"1-6 
and  for  Manchester  fine  sand: 

p  =  240  9?'  82  (1  +  9)°‘  82 

1  74  9  -1’24 

B  =  0.43  90-1-24(l  +£-) 

w0 

For  90  =  IF 

p  =  240  (1  +  9)0>  82 
B  =  0.43(1  +  9) " 1  ’  24 

where  P  is  in  psi,  B  is  in  hours  and  9  in  F  deg. 

By  substituting  computed  values  for  p  and  B  in  eq  7a,  strength  values  for 
various  temperatures  and  periods  of  stress  application  can  be  computed.  Com¬ 
puted  strength  values  plotted  on  Figures  43  and  44  are  in  good  agreement  with 
the  data  curves  shown  on  the  same  figures. 

The  increase  in  strength  of  frozen  soil  as  the  temperature  decreases  is 
generally  attributed  to  the  reduction  of  the  amount  of  unfrozen  water  in  the 
soil  and  to  increased  strength  of  the  pore  ice.  Since  the  amount  of  unfrozen 
water  in  sand  is  extremely  small  and  its  influence  on  strength  may  be  neglected, 
then  the  main  factor  affecting  change  in  strength  with  temperature  is  the  pore 
ice.  Tests  performed  by  Butkovich  (1954)  on  lake  ice  and  on  commercially 
frozen  ice  showed  a  marked  increase  in  strength  with  a  decrease  in  temperature 
which  is  in  agreement  with  the  increase  in  strength  noted  in  saturated  free  en 
sands. 


GRKK.P  OK  KROZKN  SANDS 


CONCLUSIONS 

I’lu-  primary  factors  affecting  the  strength  and  drlormation  of  a  given 
saturated  frozen  sand  under  static  load  are  temperature  ;mo  duration  of  stress. 
(Density  is  anotln  r  important  factor,  however  it  was  not  investigated  in  this  re¬ 
port.)  1  lie  result.-,  ol  this  investigation  which  are  bast'd  on  nnconfined  compres¬ 
sion  tests  performed  on  saturated  frozen  Ottawa  sand  and  Manchester  fine  sand 
show  that: 

a.  Long-term  strength  is  less  than  1  ->ir‘(>  of  instantaneous  strength.  Via- 

I.  ■  1  Oi-.in.ftl,  .  ,n  1 1 ..  1  ,  >  —  .  —  _ _  _ .  f,,«  —  ....  I.  .1. 


loc'i  strength  equation,  <ru[t  = 


for  <ruji-  =  p/log(t/B)]  ,  can  be  isi  d  to 


predict  creep  strengths  of  saturated  frozen  sand  based  on  short-term  creep 
parameters  provided  that  a  fadtor  of  safety  of  l.i  to2  is  used.  For  the  density 
and  test  conditions  the  parameters  p  and  B  can  be  considered  to  be  functions 
of  temperature  determined  approximately  for  frozen  saturated  Ottawa  sand  by: 

p  =  260  0O°-91  (1  +  f-)°‘91 

Ho 

- 1  6  0  -1.6 

B  =  0.  64  0O  1 '  (  1  +  £-) 


and  for  Manchester  fine  sar.d  by: 

0  82  0  0.  82 
p  =  24O0o°-aZ(l  +  g-) 

- 1  24  0  *  1-24 

B  =  0.  43  0O  1  4  ( 1  -  J-) 

Ho 

where  0  isinpsi,  B  is  in  hours,  9  is  Fahrenheit  degrees  below  32F  and  0Q  is  a 
reference  temperature,  >0,  F  deg. 

b.  As  the  temperature  decreases,  strength  increases.  The  long-term 
strength  of  f:  ozen  saturated  Ottawa  sand  can  be  represented  approximately  by 

ffuit  s  zo'«!  ' 

where  o-yit  is  in  psi  and  0  is  in  F  deg  f the  strength  at  32F  (0  -0)  is  very  small 
and  extremely  difficult  to  obtain;  extrapolation  is  doubtful,  so  zero  strength 
is  assumed];  and  Manchester  fine  sand  by: 

u-ult  -  40  +  40  M  3M. 

e.  Deformation  increases  with  increased  stress,  time  and  temperature. 
Unconfined  compression  strains  can  be  predicted  by: 


aa.Q  ,  ..a 

°bl  9o(7T-  +  1  ) 

L  w0 


1/K  t4, 


wh  e  ,v  e : 


applied  constant  stress 
time 
( M  -  1 )  /  M 
cr1  /w 

*F  below  32 F 

reference  value  ol  0,  F  deg 


and  the  constants  oq j  ,  K,  a  and  w  are  listed  in  Table  V. 
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(2)  Vialov's  equation 


a  t 


ll/r 


6oK(f-+  D‘ 


where  the  constants  m,  \,  u  and  k  are  listed  in  Table  IV. 


d.  For  a  constant  temperature  and  an  applied  constant  stress  the  strain 
can  be  predicted  approximately  by: 


f 


4 

4 


M- 1 

where:  ^  =  ( — j^~)  and  M  and  can  be  determined  graphically  by  a  creep  test 

lasting  only  3  hours.  Further  research  is  desirable  to  investigate  a  possible 
temperature  effect  on  M. 


e.  Strains  for  stresses  less  than  the  long-term  strength  are  composed 
of  less  than  14%  elastic  strain  (of  this  less  than  7%  is  instantaneous  and  10% 
time -dependent  elastic)  and,  further,  that  the  irreversible  strain  is  dominant 
even  at  this  stress  level. 


f.  Unconfined  compression  creep  tests  with  constant  stresses  below  the 
long-term  strength  reveal  that  the  rate  of  strain  is  directly  proportional  to  the 
reciprocal  of  time  during  stress  application  until  complete  stabilization  occurs. 

This  investigation  has  attempted  to  relate  empirically  the  effect  of  tem¬ 
perature  and  time  to  both  the  strength  and  deformation  characteristics  of  re¬ 
molded,  artificially  frozen  sai.us  subjected  to  uniaxial  compression.  The 
effect  of  density  of  the  frozen  soil,  the  soil  particle  characteristics  (i.e., 
size,  shape,  mineral  composition,  etc.),  the  water  content  (both  frozen  and 
unfrozen),  and  the  influence  of  a  more  complex  state  of  stress  are  items  that 
still  remain  for  future  experimental  investigations.  A  practical  theory  that 
will  permit  a  better  utilization  of  laboratory  and  field  test  results  in  the  de¬ 
sign  of  engineering  structures  in  frozen  soil  must  still  be  developed  before 
the  real  potential  of  frozen  soil  as  an  engineering  material  can  be  realized. 
Therefore,  the  present  investigation  can  be  viewed  only  as  an  extension  of  the 
initial  step  taken  by  ACFEL  (1952)  toward  a  practical  solution  of  engineering 
problems  in  permanently  frozen  soil. 
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Dry  unit  wt  e  =  Void  ratio 

Mass  unit  wt  P.  F.  -  Point  on  time-strain  curve  where  plastic  tlow  begins 

Vol.  of  ice  *  =  Max.  stress 

Vol.  of  roil  grains  F  =  Duration  of  test 

Percent  of  voids  filled  with  ice  V  =  Tests  performed  on  constant  stress  apparatus 
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Table  AIII.  Ice. 


Specimen 

No.  g/cm* 

Stress 
lb /in.  2 
(initial) 

%  of 
Inst. 
Strength 

Temp 
•  F 

Time 

to 

fail 

Failure 

strain 

Remarks 

Nominal  temp 

25F 

ICE- 10 

0.  914 

481* 

98.2 

25 

4.2" 

- 

LMP  failure 

ICE- 12 

0.914 

478* 

97.  6 

25 

4.  0" 

. 0000756 

ICE-22 

0.912 

510* 

104.  1 

25 

4.2" 

- 

LMP  failure 

ICE- 11 

0.915 

170 

34.  7 

25 

49  hr  Ft 

- 

Specimen  tilted 

ICE-  9 

0.914 

100 

20.  4 

25 

408  hr  Ft 

- 

ICE- 14 

0.915 

50 

10.2 

25 

791  hr  Ft 

- 

Average 

instantaneous  strength  =  490  psi 

Nominal  temp  29F 

ICE-  6 

0.  903 

259* 

96.  3 

29 

1.8" 

.‘0000585 

ICE-20 

0.915 

263* 

97.8 

29 

2.  1" 

. 000498 

ICE-21 

0.914 

285* 

105  9 

29 

2.  3" 

. 001281 

ICE-16 

0.914 

160 

59.5 

29 

2.  5  hr  Ft 

- 

ICE-15 

0.  913 

95 

35.3 

29 

627  hr  Ft 

- 

ICE-  2 

0.  914 

54 

20.  1 

29 

863  hr  Ft 

- 

ICE-  8 

0.912 

27 

10.  0 

29 

1055  hr  Ft 

- 

Average  instantaneous  strength  =  269  psi 


*  Max.  stress 
|  Duration  of  test 
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Unconfined  compressive  creep  strengths  and  strains  were  measured  for  frozen 
saturated  Ottawa  sand  (20-30)  and  Manchester  fine  sand.  The  creep  tests  were  con¬ 
ducted  at  approximate  stress  levels  of  60,  35,  20  and  5%  of  the  conventional  uncon¬ 
fined  compressive  strength.  Testing  temperatures  #ere  15,  25,  29  and  31F.  It  was 
found  that  the  unconfined  compressive  creep  strength  of  the  frozen  sand  can  be  pre¬ 
dicted  using  two  short-term,  high-stress-level  creep  tests,  using  «  ■  cj  •  (t'wili); 

<r  '  1T/,K  t  + 

that  total  strain  can  be  predicted  using  <  ■  - - - —  “j?  +  «  0  •  and 

>,0/00)  +  1]  J  * 

that  for  stresses  below  the  long-term  strength,  the  strain  rate  is  directly  proportion¬ 
al  to  the  reciprocal  of  time  during  stress  action  until  complete  stabilization  occurs. 

(*  a  strain  rate  1  hour  after  stress  is  applied;  t  =  time;  -  (M-l)/M,  where 
M  x  g.l/vi  and  w  is  a  constant  for  each  material;  <%  *  stress  at  0  =  0;  0  »  tempera. ui 
in  degrees  below  freezing  point  of  water;  Oq  =  a  constant  reference  value  o.  0,  u  and 
K  are xonstants;  «  •*  initial  instantaneous  strain.) 
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